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observed B+ and B° decays, B{B+ -> = (0.59 ± 0.15 ± 0.04) x 10~ 6 , and B(B° -> 
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The decay of a & quark into an s quark and two muons 
(b — > sfi + n~) is a flavor-changing neutral current pro- 
cess, forbidden at tree level in the standard model (SM) 
but allowed through highly suppressed internal loops. 
New physics could manifest itself in a larger branching 
fraction, a modified dimuon mass distribution, or an- 
gular distributions of the decay products different from 
that predicted by the SM 0, [2, d, In this pa- 
per, we report branching ratio measurements of exclu- 
sive decays, where the s quark hadronizes into a single 
meson. Reconstructing a final state with only three or 
four charged final state particles results in smaller back- 
grounds than those expected in a search for inclusive 
B — > X s ^l + n~ decays. The rare decays, B + — » [i + n~K + 
and B° — > n + \i~ K* (892)° have been observed at the B 
factories p||f|, with branching fractions of O(10~ 6 ). con- 
sistent with SM predictions @, B i, (M El El, El El- 
The analogous decay in the B° system, B° — > [i + [i~(f> 7 
has a predicted branching ratio of 1.6 x 10 -6 [l5j|, but has 
not yet been observed [H E3 • We report branching ratio 
measurements from 924 pb -1 of Collider Detector at Fer- 
milab (CDF) Run II data that pave the way for future 
studies of larger datasets from which kinematic distri- 
butions may be measured. We search for B — > [i + [i~h 
decays, where B stands for B + , B°, or B®, and h stands 
for K+, K* (892)°, or </>, respectively. The if* (892)°, re- 
ferred to as K*° throughout this paper, is reconstructed 
in the K* a — > K + ir~, decay mode, and the 4> meson 
is reconstructed as (f> — > K + K~ . We measure branch- 
ing ratios relative to B — > J/iph decays, followed by 
J/ip — > decays, resulting in the same final state 

particles as the rare decay modes. Many systematic un- 
certainties cancel in the relative branching ratios: 



B(B^n+n-h) _ N, 
B(B -» J/iph) ~ N 



fi+fj.-h e J/i>h 



where N^+u- h is the observed number of B — ► h de- 

cays, Nj/^ h is the observed number of B — > J/iph decays, 
while ej/^h and e^+^-h ar e the combined selection effi- 
ciency and acceptance of the experiment for B — > J/iph 
and _B — ► fj, + fi~h respectively. Throughout this report, 
charge conjugate modes are implicitly included. 

CDF II is a general purpose detector, located at the 
Tevatron pp collider jla |. Charged particle trajectories 
(tracks) are detected by the tracking system comprised of 
a seven-layer double-sided silicon microstrip detector and 
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a drift chamber, both in a 1.4 T axial magnetic field. The 
silicon detector [n| ranges in radius from 1.3 to 28 cm, 
and has a single- hit resolution of approximately 15 fim. 
The drift chamber [2(| provides up to 96 measurements 
from radii of 40 to 137 cm with a single-hit resolution of 
approximately 180 /im. A muon chamber identification 
system of plastic scintillators and drift chambers [2l| is 
located on the exterior of the detector with a central 
part covering \tj\ < 0.6, and an extended region cover- 
ing 0.6 < | r) | < 1.0, where r\ is the pseudorapidity [22I ]. 
In the central(extended) region, muons are detected if 
their transverse momentum component, pt, is greater 
than 1.5(2.0) GeV/c. Events are selected with a three- 
level trigger system. The first trigger level requires the 
presence of two charged particles with pr > 1.5GeV/c 
(hi < 0.6) or p T > 2.0GeV/c (0.6 < \r]\ < 1.0), matched 
to track segments in the muon chambers to form muon 
candidates. At the second level, a more restrictive se- 
lection is made by requiring that the muon candidates 
have opposite charge, and that their opening angle in the 
plane transverse to the beamline is less than 120°. At the 
third trigger level, the event is fully reconstructed. The 
trajectories of the muon candidates in the silicon detec- 
tor are required to intersect at a point which is displaced 
transversely from the beamline by at least 100 /im. 

The offline selection begins with a suppression of ran- 
dom combinations of tracks satisfying the selection re- 
quirements (combinatoric background), by requiring that 
all tracks have px > 0.4GeV/c and match hits from at 
least three layers of the silicon detector. The trajecto- 
ries of a pair of muon candidates that satisfy the trigger 
requirements, and the tracks that form the hadron can- 
didate, are fitted with the constraint that they originate 
from a single vertex in 3-dimensional space to form a 
B candidate. The \ 2 probability of the fit is required 
to be greater than 10~ 3 . A K + candidate is a track 
assigned the charged kaon mass, and the K*° (<jf>) candi- 
dates are formed from oppositely charged pairs of tracks 
whose invariant mass must lie within 50 (10)MeV/c 2 of 
the K*° ((f)) mass. For all particles, we use the world av- 
erage values tabulated in Ref. HI . The ambiguity of the 
mass assignment in B° — > fi + fi~K*° decays is handled 
by choosing the combination whose K + ir~ mass is closer 
to the K*° mass. In reconstructing the B candidates, the 
meson containing a strange quark, h, is required to have 
Pr(h) > l.OGeV/c, and the B candidate is required to 
have pr(B) > 4.0GeV/c. We require that the distance 
of closest approach between the flight path of the B can- 
didate and the beamline, |do(i?)|, is less than 120 /im, 
to reduce the combinatoric background with no loss of 
signal. 

Normalization mode candidates are identified by hav- 
ing a dimuon invariant mass within 50MeV/c 2 of the 
J/ip mass, yielding approximately 12 000 B + — » J /ipK + , 
4400 B° -> J/i;K* Q , and 800 B° s -> J/ifxf) decays. The 
B candidate mass distributions of all modes are com- 
patible with a Gaussian of width a = 20MeV/c 2 . To 
reduce backgrounds from B decays to mesons contain- 
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ing c quarks, several vetoes are applied to B — > 
candidates, listed below. We eliminate candidates with a 
dimuon mass near the J/ipW : 2.9 < m w < 3.2GeV/c 2 
or 3.6 < ra w < 3.75GeV/c 2 respectively. B -> J/ip^h 
decays followed by the radiative decay of the J/ip^ into 
two muons and a photon that is not reconstructed, may 
have a dimuon mass that passes the J/ip^'' veto. We re- 
ject these events making use of the correlation between 
the candidate's invariant mass and the dimuon invariant 
mass, and require {(m^h ~ rn B ) - (m^ - rrij/^)\ > 
100MeV/c 2 . B -> J/^h decays with one hadron 
misidentified as a muon form a potential background to 
the rare decay search. We reject this class of background 
by requiring that all combinations of tracks of the can- 
didate have an invariant mass that differs by at least 
40MeV/c 2 from the J/ip^ mass. We also reject candi- 
dates with track pairs with a mass within ±25MeV/c 2 
of the D° — > K~ir + decay, or track triplets compatible 
with a mass within ±25MeV/c 2 of the D+ -> A^7r + 7r+ 
or Z?+ — > K + K~tt + decays. 

We further improve the signal to background ratio by 
an optimization of the selection based on discriminating 
variables. For this purpose, Monte Carlo (MC) simula- 
tions of the signal processes are used. We generate sin- 
gle b hadrons using the transverse momentum spectrum 
from B — > J/ipX, measured by CDF (1811. Decays of all b 
hadrons are simulated using evtgen [24[, with lifetimes 
from Ref. [23j, allowing for a decay width difference be- 
tween the B° s mass eigenstates, AT /V = 0.12 ± 0.06 [25[. 
The decays B° — > J/ipK*° and B° — ► J/ip<j) are simu- 
lated according to the polarization amplitudes measured 
by CDF . The dynamics of the rare decay pro- 
cesses are simulated according to the calculations of Ali 
et at. 11]. For B° — > fi + fi~cj), we assume a mixture of 
50% CP-even and 50% CP-odd states. The interactions 
of final state particles are simulated using a GEANT [27| 
model of the CDF II detector, digitized into the CDF 
event format, and reconstructed using the same software 
as in the processing of collision data. The detector sim- 
ulation includes a full emulation of the CDF trigger sys- 
tem. We find that the transverse momentum spectrum of 
the B mesons of this measurement are on average higher 
than in the simulation. This can be explained by the 
presence of a B baryon component with low transverse 
momentum in the measurement used as input spectrum 
of the simulation [18J. To correct for this, the simulated 
events are weighted by a polynomial function, which is 
obtained from the ratio of data to MC transverse mo- 
mentum distribution of the B meson. 

We tighten the candidate selection to obtain the small- 
est expected statistical uncertainty on the measurement 
of the B — > n + ii~h branching ratios by finding the selec- 
tion with the largest value of N sig / ^jN S i g + Nbkg ■ The 
expected number of signal events, N S ig, is determined 
by scaling the MC yields of the rare decay modes to 
the yields of the normalization channels in the data, 
and correcting by the relative branching ratios and se- 
lection efficiencies. This procedure requires estimates of 



the branching ratios of the rare decay modes. For the 
B + — > p, + fi~K + and B° — » fi + ^K*° decays we use 
the world average measured branching ratios [23| . while 
for B® — > /i + /i~(/> we use the theoretical estimate (l5j . 
The expected number of background events, N^kg, is es- 
timated by extrapolating the candidate yield with an 
invariant mass (180 — 300) MeV/c 2 higher than the B 
mass, to the signal region. Candidates with an invari- 
ant mass lower than the B mass are not suitable for 
background estimates, since they contain partially re- 
constructed B decays, not expected in the signal region. 
We find good discriminating power between signal and 
background from the following three quantities: t/at, a, 
and I. The significance of the proper decay time, t/at, 
is defined as the proper decay time of the B candidate, 
divided by its uncertainty. The angle a is defined as 
the difference in angle between the B candidate's mo- 
mentum vector and the vector from the primary vertex 
to the \i[ih vertex. The isolation, /, is defined as the 
transverse momentum carried by the B meson candidate 
divided by the transverse momentum of all tracks in a 
cone of AR — \J Arf + A(p 2 = 1.0 around the direc- 
tion of the B meson candidate, including those of the 
B candidate itself. Here Ar] is the difference in pseudo- 
rapidity of the B candidate and each track, and A</> is 
the difference in their azimuthal angles. Scanning dif- 
ferent combinations of selection thresholds, we find that 
the optimal values are very similar for the three rare 
decay modes: t/at > 14, a < 60mrad, and I > 0.6. 
Applying the optimized selection requirements to the 
normalization channels yields the distributions shown in 
Fig. [TJ The shape of the combinatoric background is 
estimated from a sample with poor vertex quality (x 2 
probability < 10~ 3 ). The invariant mass distribution of 
this poor vertex quality sample is fitted with a Gaus- 
sian distribution for the remaining signal contribution 
and an exponential plus a constant to model the back- 
ground. We estimate the background in the signal region 
of the tight selection using this functional form, normal- 
ized to the number of candidates with invariant mass 
(60 - 180) MeV/c 2 higher than the B mass. Using this 
method, we find 6361 ± 82 B+ -> J/t/>A+, 2423 ± 52 
B° -> J/tpK*°, and 431 ± 22 B° s -> J/^<\> decays. The 
invariant mass distributions for the rare decay modes are 
shown in Fig. O In the ±40MeV/c 2 window around the 
B mass we find 90 B+ -> n+fj,~K + , 35 B° -> fj+^-K* , 
and 11 B® — > /i + 4> candidates. From simulation stud- 
ies we verified that the present selection accepts decays 
for every kincmatically possible dimuon invariant mass, 
with an efficiency difference not exceeding a factor two, 
except for the windows around the J/^ 1 ' that have been 
explicitly excluded. 

For the signal modes, we evaluate three sources of 
background: the combinatoric background, hadrons 
misidentified as muons, and hadrons with misassigned 
mass. As in the normalization mode, the shape of the 
combinatoric background is estimated from a sample 
with poor vertex quality (x 2 probability < 10 -3 ). The 
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Invariant mass spectra of (a) J/ipK + , (b) J/ipK*° and (c) J/i/xj) candidates after applying optimized selection require- 



invariant mass distribution of these candidates is fitted to 
an exponential plus a constant. This functional form is 
used to extrapolate the number of candidates with invari- 
ant mass (60 — 180) MeV/c 2 higher than the B mass to 
the signal region. The resulting backgrounds, and their 
statistical uncertainties, consist of 44.3 ± 5.8, 16.3 ± 3.6, 
and 3.1 ± 1.5 candidates for the B + — > /i + /i~A' + , 
B° — > fi + [i~K* , and B° — > fi + \i~ <f) decay modes, re- 
spectively. Backgrounds from hadrons misidentified as 
muons are estimated from simulation. Final state parti- 
cles from simulated charmless B decays are weighted with 
the muon misidentification probabilities. The misiden- 
tification probabilities for charged kaons and pions are 
measured from a sample of D* + — > D°ir + followed by 
D° — > K~n + decays. We estimate a net background 
of 1.0 events for the B + — * (i + [j,~K + channel, primar- 
ily coming from misidentified B + — > _ftT + 7r~7r + decays. 
For the £?° and B® modes, this background is negligible. 
From simulation we determine the background contribu- 
tion that arises from assigning the wrong masses to the 
hadrons. We estimate 0.4 B° — > ^, + fi~K*° decays which 
are reconstructed as B° s — > and 0.2 — > /j + /i~</> 

decays reconstructed as S — > fi~ K*°. 

We calculate the signal yield by subtracting the pre- 
dicted background from the number of candidates in the 
signal window. We find an excess in the B signal region 
in all three channels, and determine the significance by 
calculating the Poisson probability for the background 
to fluctuate to the number of observed events or higher, 
taking into account the uncertainty on the background. 
We find an equivalent Gaussian significance of 4.5, 2.9, 
and 2.4 standard deviations, respectively for the B + , B , 
and Bg modes. 

Table U lists the systematic uncertainties associated 
with the relative selection efficiency, discussed in fur- 
ther detail below. To take into account the uncertainty 
on the relative efficiency due to uncertainties in the dy- 
namics of the rare decays, we apply the weak form fac- 
tors from Ref. [28|, and evaluate the largest differ- 
ence between the reconstruction efficiency and the cen- 



tral value. We find an uncertainty of 3.1% or less. The 
uncertainty related to the px(B) spectrum is evaluated 
from the change in relative efficiency when the three pa- 
rameters in the pt{B) weighting function are varied by 
one standard deviation of their values determined from 
fits to data, taking into account correlations. The dif- 
ference reaches 1.4% in the B® mode. The muon trigger 
efficiency close to the 1.5GeV/c pr threshold is poorly 
known. We find a change in the relative efficiency reach- 
ing 1.3% if the minimum pr is varied by ±100MeV/c, 
covering the range over which the trigger efficiency rises 
from zero to one. The final state particles of the rare 
decay modes have approximately 10% more low momen- 
tum tracks than those of the normalization channels. The 
simulation models the track reconstruction efficiency to 
an accuracy of 2% in the pr range of 0.4 — 1.5 GeV/c. We 
therefore assign a systematic uncertainty on the relative 
efficiency of 2%xl0%=0.2%. We find an uncertainty of 
8.3% on the relative efficiency of B® decays due to the 
unknown fraction of the short-lived CP-even state in the 
rare decay mode. We assume that the CP-even frac- 
tion is 0.5, and assess systematic uncertainties for the 
extremes of the fraction valued at 1.0 and 0.0. In addi- 
tion, the uncertainty on Ar/r contributes another 2.6%, 
resulting in a total uncertainty of 8.7% associated with 
the B® decay width difference. This is the largest sys- 
tematic uncertainty on the relative efficiency of the B® 
mode. We evaluate the effect of the uncertainty on the 
fraction of J/tp mesons produced with a longitudinal po- 
larization by varying the fraction measured at CDF [26j 
by ±lcr. The effect is 0.6% for B° -► J/^K*° and 0.1% 
for B° s -> J/ip<p. 

The statistical uncertainties of yields in the normal- 
ization channels arc included as systematic uncertainties, 
which range from 1.3% for the B + channel to 5.1% for the 
B° channel. Cabibbo-suppressed B + — > J/^tt + decays 
contribute 0.1% to the yield of the normalization channel 
for the B + channel. We introduce the full effect into our 
estimate of the systematic uncertainty without correct- 
ing the result. The relative efficiencies between signal 
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TABLE I: Systematic uncertainties on the relative efficiency 
quoted in percent. 



Channel B+ B u BT 

Theory model 1.5 3.1 1.6 

p T {B) spectrum 0.6 1.3 1.4 

Trigger turn-on 1.3 1.3 1.2 

Low momentum hadrons 0.2 0.2 0.2 

Bg decay width difference — — 8.7 

Polarization — 0.6 0.1 

Norm, channel statistics 1.3 2.1 5.1 

B + — > J/ip-K + contribution 0.1 — — 

MC statistics 1.6 2.6 2.2 

Tbtal 2.9 5.0 10.6 



TABLE If: Summary of systematic uncertainties quoted in 
percent. 



Channel 


B+ 


B u 


B s 


Total rel. eff. uncertainty 


2.9 


5.0 


10.6 


B(J/i>^IM + fT) 


1.0 


1.0 


1.0 


Background prediction 


5.2 


3.1 


10.2 


B(B -> J/ij>h) 


3.5 


4.5 


35.5 



and normalization channels, of 0.71 ± 0.01, 0.74 ± 0.02, 
and 0.84 ± 0.02 for the B+ , B°, and B° decays, respec- 
tively, have uncertainties reaching 2.6%, that arise from 
the finite size of the MC samples. 

The predicted background values depend on the shape 
of the background function. We evaluate the change in 
the background yield calculation when using a sample 
that is similar to the one resulting from the optimal se- 
lection, but has less stringent requirements on the three 
optimization variables instead of the default sample with 
poor vertex quality. We compare two functional models 
for the background description, the default exponential 
plus a constant, and a simpler linear extrapolation func- 
tion. We find a systematic uncertainty on the background 
prediction reaching 10% for the B® channel. 

To calculate the systematic uncertainty on the rela- 
tive branching fraction, the relative efficiency uncertainty 
is summed in quadrature with the uncertainty of the 
J/tp — » branching ratio, and with the systematic 

uncertainty on the number of signal candidates. 

To calculate the absolute branching fractions, we use 
the world average branching fractions of the normaliza- 
tion channels [23|]. These branching fractions have un- 
certainties of 3.5%, 4.5%, and 35.5%, respectively, for 
B + , B°, and B°, which are added in quadrature to the 
systematic uncertainties on the relative branching ratios. 
A summary of the systematic uncertainties is given in 
Table HI 

We calculate relative branching ratios using Eq. fl}, 
and find B(B+ -> /z+ fi~ K + )/B(B + -> J/ipK + ) = 
(0.59±0.15±0.03) x 10~ 3 , B{B° -> fi+ >- K* )/ 'B(B° -» 



J/ifiK* ) = (0.61 ± 0.23 ± 0.07) x lO" 3 , and B{B° S -» 
fi+[i-<f>)/B(B° -» J/ip<f>) = (1.23 ±0.60 ±0.14) x 10~ 3 , 
where the first uncertainty is statistical, and the sec- 
ond is systematic. We use the world-average branch- 
ing ratios of the B° and B + normalization channels [23| , 
resulting in the following absolute branching fractions: 
B(B+ -> fi+fi-K+) = (0.59 ±0.15 ±0.04) x 10~ 6 , and 
B(B° -> n+fi-K* ) = (0.81 ± 0.30 ± 0.10) x 10~ 6 . To 
obtain an absolute branching ratio for the B® — > ^ \jT '<f> 
decay, we use B{B° S -> J/^) = (1.38 ± 0.49) x 10" 3 , 
obtained from correcting the CDF measurement [3(3] for 
the current value of f s /fd ,23], the B° s to B° produc- 
tion ratio, resulting in B{B° S -» n + V~<f>) = (1.70 ± 0.82 ± 
0.64) x 10 -6 . We find a significance of only 2.4 standard 
deviations for the B® fi + [i~(j) decay mode. There- 
fore, we choose to set a limit on this decay. We use a 
Bayesian integration assuming a flat prior [3l| . and find 
B(B° S -> ^+n-(f))/B{B Q s -> J/ip<f>) < 2.6(2.3) x 10~ 3 at 
the 95(90)% confidence level (C.L.). We also set an upper 
limit on the absolute branching ratio, taking into account 
the uncertainty of the branching ratio of the normaliza- 
tion channel [321 ] . We assume that the prior probability 
density function representing the uncertainty in the nor- 
malization channel is a log normal distribution with a 
mean equal to the central value, and a width parame- 
ter equal to the quoted uncertainty of 35.5%. We obtain 
B(B° -> /i+m~0) < 6.0(5.0) x 10~ 6 at 95(90)% C.L. The 
main ingredients and the results are summarized in Ta- 
ble Ell 

In conclusion, our measurements of the B + and B° rare 
decay modes are consistent with the SM predictions, and 
with previous measurements 0, @ . The relative limit on 
£>!? — > is consistent with the SM predictions, and 

is the most stringent to date. 
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FIG. 2: Invariant mass spectra of (a) fi~ K + , (b) fi + fi~ K*° and (c) fi + fj,~(j> candidates. The superimposed curves are a sum 
of a single Gaussian with a width of 20 MeV/c 2 representing the signal, and a curve representing the background as determined 
by the procedure described in the text. The curves are not drawn for masses below the signal window, since they are not 
expected to predict the background level where partially reconstructed B decays contribute. 



TABLE III: Summary of the main ingredients and results of this analysis. Where one uncertainty is quoted, the uncertainty is 

of statistical nature. Where two uncertainties are quoted, the first is statistical, and the second is systematic. 

T3ecaymode B+ -> /i + ^ K+ B u -> y+ fi~ A'* u B s ' -> n + fi~ (j> 

N oba 90 35 11 

N bkg 45.3 ±5.8 16.5 ±3.6 3.5 ±1.5 

Nsig 44.7 ±5.8 18.5 ±3.6 7.5 ±1.5 

Gaussian significance 4.5 a 2.9 a 2.4 a 

N J/i>h 6361 ±82 2423 ± 52 431 ± 22 

S+M-h/ e J/*h 0.71 ±0.01 0.74 ±0.02 0.84 ±0.02 

Rel B x 10 3 0.59 ± 0.15 ± 0.03 0.61 ± 0.23 ± 0.07 1.23 ± 0.60 ± 0.14 

Abs B x 10 6 0.59 ± 0.15 ± 0.04 0.81 ± 0.30 ± 0.10 1.70 ± 0.82 ± 0.64 

Rel B 95(90)%C.L. limit xlO 3 2.6(2.3) 

Abs B 95(90)%C.L. limit xlO 6 6.0(5.0) 



[1] P. Colangelo, F. De Fazio, R. Ferrandes, and T. N. Pham, 

Phys. Rev. D 73, 115006 (2006). 
[2] F. Kruger and J. Matias, Phys. Rev. D 71, 094009 (2005). 
[3] T. M. Aliev, A. Ozpineci, and M. Savci, Eur. Phys. J. C 

29, 265 (2003). 
[4] D. A. Demir, K. A. Olive, and M. B. Voloshin, Phys. 

Rev. D 66, 034015 (2002). 
[5] B. Aubert et al. (BABAR Collaboration), Phys. Rev. D 

73, 092001 (2006). 
[6] A. Ishikawa et al. (Belle Collaboration), Phys. Rev. Lett. 

91, 261601 (2003). 
[7] G. Burdman, Phys. Rev. D 52, 6400 (1995) 
[8] P. Colangelo, F. De Fazio, P. Santorelli, and E. Scrimieri, 

Phys. Rev. D 53, 3672 (1996) [Erratum-ibid. D 57, 3186 

(1998)] 

[9] T. M. Aliev, A. Ozpineci, and M. Savci, Phys. Rev. D 
56, 4260 (1997) 
[10] D. Melikhov, N. Nikitin, and S. Simula, Phys. Rev. D 57, 
6814 (1998) 

[11] A. Ali, P. Ball, L. T. Handoko, and G. Hiller, Phys. Rev. 
D 61, 074024 (2000). 



[12] H. M. Choi, C. R. Ji, and L. S. Kisslinger, Phys. Rev. D 

65, 074032 (2002) 
[13] A. Ali, E. Lunghi, C. Greub, and G. Hiller, Phys. Rev. 

D 66, 034002 (2002) 
[14] C. H. Chen and C. Q. Ceng, Phys. Rev. D 66, 094018 

(2002) 

[15] C. Q. Ceng and C. C. Liu, J. Phys. G 29, 1103 (2003). 
[16] V. M. Abazov et al. (DO Collaboration), Phys. Rev. D 

74, 031107 (2006). 
[17] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 65, 

111101 (2002). 

[18] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71, 
032001 (2005). 

[19] A. Sill et al, Nucl. Instrum. Methods A 447, 1 (2000). 
[20] A. A. Affolder et al, Nucl. Instrum. Methods A 526, 249 
(2004). 

[21] G. Ascoli et al, Nucl. Instrum. Methods A 268, 33 
(1988). 

[22] The CDF reference frame uses cylindrical coordinates, 
where and <f> are the polar and azimuthal angles with 
respect to the proton beam. Pseudorapidity rj is defined 



9 



as -ln(tan(0/2)). 
[23] W.-M. Yao et al., J. Phys. G 33, 1 (2006). 
[24] D. J. Lange, Nucl. Instrum. Methods A 462, 152 (2001). 
[25] I. Dunietz, R. Fleischer, and U. Nierste, Phys. Rev. D 

63, 114015 (2001). 
[26] D. Acosta et al. (CDF Collaboration), Phys. Rev. Lett. 

94, 101803 (2005). 
[27] R. Brun et al., CERN Program Library Long Writeup 

W5013 (1994). 

[28] D. Melikhov and B. Stech, Phys. Rev. D 62, 014006 
(2000). 



[29] P. Colangelo, F. De Fazio, and P. Santorelli, Phys. Rev. 

D 51, 2237 (1995). 
[30] F. Abe et al, (CDF Collaboration), Phys. Rev. D 54, 

6596 (1996). 

[31] J. Heinrich, C. Blocker, J. Conway, L. Demor- 
tier, L. Lyons, G. Punzi, and P. K. Sinervo, 
|arXiv:physics/0409l"29j 

[32] Note that the uncertainty in the normalization mode's 
branching ratio was not taken into account in fit} ]. 



Search for the Rare Decays B+ -> fi + fi~K + , B° -> / u + / u- J ftT*(892) , and B° -> at 

CDF 

T. Aaltonen, 24 J. Adelman, 14 T. Akimoto, 56 M.G. Albrow, 18 B. Alvarez Gonzalez, 12 S. Amerio", 44 D. Amidei, 35 
A. Anastassov, 39 A. Annovi, 20 J. Antos, 15 G. Apollinari, 18 A. Apresyan, 49 T. Arisawa, 58 A. Artikov, 16 
W. Ashmanskas, 18 A. Attal, 4 A. Aurisano, 54 F. Azfar, 43 P. Azzurri 8 , 47 W. Badgett, 18 A. Barbaro-Galtieri, 29 
V.E. Barnes, 49 B.A. Barnett, 26 V. Bartsch, 31 G. Bauer, 33 P.-H. Beauchemin, 34 F. Bedeschi, 47 P. Bednar, 15 
D. Bcecher, 31 S. Behari, 26 G. Bellettini 9 47 J. Bellinger, 60 D. Benjamin, 17 A. Beretvas, 18 J. Bcringer, 29 A. Bhatti, 51 

M. Binkley, 18 D. Biscllo", 44 I. Bizjak, 31 R.E. Blair, 2 C. Blocker, 7 B. Blumenfeld, 26 A. Bocci, 17 A. Bodek, 50 
V. Boisvcrt, 50 G. Bolla, 49 D. Bortoletto, 49 J. Boudreau, 48 A. Boveia, 11 B. Brau, 11 A. Bridgeman, 25 L. Brigliadori, 44 
C. Brombcrg, 36 E. Brubaker, 14 J. Budagov, 16 H.S. Budd, 50 S. Budd, 25 K. Burkett, 18 G. Busetto", 44 P. Bussey 1 , 22 
A. Buzatu, 34 K. L. Byrum, 2 S. Cabrera? 5 , 17 C. Calancha, 32 M. Campanelli, 36 M. Campbell, 35 F. Canelli, 18 

A. Canepa, 46 D. Carlsmith, 60 R. Carosi, 47 S. CarrihV, 19 S. Carron, 34 B. Casal, 12 M. Casarsa, 18 A. Castro*, 6 
P. Catastini r , 47 D. Cauz™, 55 V. Cavaliere r , 47 M. Cavalli-Sforza, 4 A. Cerri, 29 L. Cerrito", 31 S.H. Chang, 28 

Y.C. Chen, 1 M. Chcrtok, 8 G. Chiarclli, 47 G. Chlachidzc, 18 F. Chlebana, 18 K. Cho, 28 D. Chokheli, 16 J.P. Chou, 23 

G. Choudalakis, 33 S.H. Chuang, 53 K. Chung, 13 W.H. Chung, 60 Y.S. Chung, 50 C.I. Ciobanu, 45 M.A. Ciocci r , 47 
A. Clark, 21 D. Clark, 7 G. Compostella, 44 M.E. Convery, 18 J. Conway, 8 K. Copic, 35 M. Cordelli, 20 G. Cortiana", 44 
D.J. Cox, 8 F. Crescioli*, 47 C. Cuenca AlmenarP, 8 J. Cuevas m , 12 R. Culbertson, 18 J.C. Cully, 35 D. Dagenhart, 18 
M. Datta, 18 T. Davies, 22 P. de Barbara, 50 S. De Cecco, 52 A. Deisher, 29 G. De Lorenzo, 4 M. Dcll'Orso 9 , 47 
C. Deluca, 4 L. Demortier, 51 J. Deng, 17 M. Deninno, 6 P.F. Derwent, 18 G.P. di Giovanni, 45 C. Dionisi", 52 

B. Di Ruzza™, 55 J.R. Dittmann, 5 M. D'Onofrio, 4 S. Donate, 47 P. Dong, 9 J. Donini, 44 T. Dorigo, 44 S. Dube, 53 
J. Efron, 40 A. Elagin, 54 R. Erbacher, 8 D. Errcde, 25 S. Errede, 25 R. Eusebi, 18 H.C. Fang, 29 S. Farrington, 43 

W.T. Fedorko, 14 R.G. Feild, 61 M. Feindt, 27 J.P. Fernandez, 32 C. Ferrazza 8 , 47 R. Field, 19 G. Flanagan, 49 R. Forrest, 8 
M. Franklin, 23 J.C. Freeman, 18 I. Furic, 19 M. Gallinaro, 52 J. Galyardt, 13 F. Garberson, 11 J.E. Garcia, 47 
A.F. Garfinkel, 49 K. Genser, 18 H. Gerberich, 25 D. Gerdes, 35 A. Gessler, 27 S. Giagu 11 , 52 V. Giakoumopoulou, 3 
P. Giannetti, 47 K. Gibson, 48 J.L. Gimmcll, 50 CM. Ginsburg, 18 N. Giokaris, 3 M. Giordani 1 ", 55 P. Giromini, 20 
M. Giunta" 3 , 47 G. Giurgiu, 26 V. Glagolev, 16 D. Glcnzinski, 18 M. Gold, 38 N. Goldschmidt, 19 A. Golossanov, 18 
G. Gomez, 12 G. Gomez-Ceballos, 33 M. Goncharov, 54 O. Gonzalez, 32 I. Gorelov, 38 A.T. Goshaw, 17 K. Goulianos, 51 

A. Gresele", 44 S. Grinstcin, 23 C. Grosso-Pilcher, 14 R.C. Group, 18 U. Grundlcr, 25 J. Guimaraes da Costa, 23 
Z. Gunay-Unalan, 36 C. Haber, 29 K. Hahn, 33 S.R. Hahn, 18 E. Halkiadakis, 53 B.-Y. Han, 50 JY. Han, 50 R. Handler, 60 

F. Happacher, 20 K. Hara, 56 D. Hare, 53 M. Hare, 57 S. Harper, 43 R.F. Harr, 59 R.M. Harris, 18 M. Hartz, 48 
K. Hatakeyama, 51 J. Hauser, 9 C. Hays, 43 M. Heck, 27 A. Heijboer, 46 B. Heinemann, 29 J. Heinrich, 46 C. Henderson, 33 
M. Hcrndon, 60 J. Heuser, 27 S. Hewamanage, 5 D. Hidas, 17 C.S. Hill , 11 D. Hirschbuehl, 27 A. Hockcr, 18 S. Hou, 1 
M. Houlden, 30 S.-C. Hsu, 10 B.T. Huffman, 43 R.E. Hughes, 40 U. Husemann, 61 J. Huston, 36 J. Incandela, 11 

G. Introzzi, 47 M. Iori", 52 A. Ivanov, 8 E. James, 18 B. Jayatilaka, 17 E.J. Jcon, 28 M.K. Jha, 6 S. Jindariani, 18 

W. Johnson, 8 M. Jones, 49 K.K. Joo, 28 SY. Jun, 13 J.E. Jung, 28 T.R. Junk, 18 T. Kamon, 54 D. Kar, 19 P.E. Karchin, 59 
Y. Kato, 42 R. Kcphart, 18 J. Kcung, 46 V. Khotilovich, 54 B. Kilminster, 40 D.H. Kim, 28 H.S. Kim, 28 J.E. Kim, 28 
M.J. Kim, 20 S.B. Kim, 28 S.H. Kim, 56 Y.K. Kim, 14 N. Kimura, 56 L. Kirsch, 7 S. Klimcnko, 19 B. Knuteson, 33 

B. R. Ko, 17 S.A. Koay, 11 K. Kondo, 58 D.J. Kong, 28 J. Konigsberg, 19 A. Korytov, 19 A.V. Kotwal, 17 M. Kreps, 27 
J. Kroll, 46 N. Krumnack, 5 M. Kruse, 17 V. Krutelyov, 11 T. Kubo, 56 T. Kuhr, 27 N.P. Kulkarni, 59 M. Kurata, 56 

Y. Kusakabe, 58 S. Kwang, 14 A.T. Laasanen, 49 S. Lami, 47 S. Lammel, 18 M. Lancaster, 31 R.L. Lander, 8 K. Lannon, 40 
A. Lath, 53 G. Latino r , 47 I. Lazzizzera", 44 T. LeCompte, 2 E. Lee, 54 S.W. Lee°, 54 S. Leone, 47 S. Levy, 14 J.D. Lewis, 18 
C.S. Lin, 29 J. Linacre, 43 M. Lindgren, 18 E. Lipeles, 10 A. Lister, 8 D.O. Litvintsev, 18 C. Liu, 48 T. Liu, 18 
N.S. Lockyer, 46 A. Loginov, 61 M. Loreti", 44 L. Lovas, 15 R.-S. Lu, 1 D. Lucchesi", 44 J. Lueck, 27 C. Luci", 52 
P. Lujan, 29 P. Lukcns, 18 G. Lungu, 51 L. Lyons, 43 J. Lys, 29 R. Lysak, 15 E. Lytken, 49 P. Mack, 27 D. MacQueen, 34 
R. Madrak, 18 K. Maeshima, 18 K. Makhoul, 33 T. Maki, 24 P. Maksimovic, 26 S. Malde, 43 S. Malik, 31 G. Manca, 30 
A. Manousakis-Katsikakis, 3 F. Margaroli, 49 C. Marino, 27 CP. Marino, 25 A. Martin, 61 V. Martin 1 , 22 M. Martinez, 4 
R. Martinez-Ballarin, 32 T. Maruyama, 56 P. Mastrandrea, 52 T. Masubuchi, 56 M.E. Mattson, 59 P. Mazzanti, 6 
K.S. McFarland, 50 P. Mclntyre, 54 R. McNulty' 1 , 30 A. Mehta, 30 P. Mehtala, 24 A. Mcnzione, 47 P. Merkel, 49 

C. Mesropian, 51 T. Miao, 18 N. Miladinovic, 7 R. Miller, 36 C. Mills, 23 M. Milnik, 27 A. Mitra, 1 G. Mitselmakher, 19 

H. Miyake, 56 N. Moggi, 6 C.S. Moon, 28 R. Moore, 18 M.J. Morello*, 47 J. Morlok, 27 P. Movilla Fernandez, 18 
J. Miilmenstadt, 29 A. Mukherjee, 18 Th. Muller, 27 R. Mumford, 26 P. Murat, 18 M. Mussini*, 6 J. Nachtman, 18 
Y. Nagai, 56 A. Nagano, 56 J. Naganoma, 58 K. Nakamura, 56 I. Nakano, 41 A. Napier, 57 V. Necula, 17 C. Neu, 46 



M.S. Neubauer, 25 J. Nielsen 6 , 29 L. Nodulman, 2 M. Norman, 10 O. Norniella, 25 E. Nurse, 31 L. Oakes, 43 S.H. Oh, 17 
Y.D. Oh, 28 I. Oksuzian, 19 T. Okusawa, 42 R. Oldcman, 30 R. Orava, 24 K. Osterberg, 24 S. Pagan Griso", 44 
C. Pagliarone, 47 E. Palencia, 18 V. Papadimitriou, 18 A. Papaikonomou, 27 A. A. Paramonov, 14 B. Parks, 40 
S. Pashapour, 34 J. Patrick, 18 G. Pauletta™, 55 M. Paulini, 13 C. Pans, 33 D.E. Pellett, 8 A. Penzo, 55 T.J. Phillips, 17 
G. Piaccntino, 47 E. Pianori, 46 L. Pinera, 19 K. Pitts, 25 C. Plager, 9 L. Pondrom, 60 O. Poukhov, 16 N. Pounder, 43 

F. Prakoshyn, 16 A. Pronko, 18 J. Proudfoot, 2 F. PtohosV 8 E. Pueschel, 13 G. Punzi«, 47 J. Pursley, 60 
J. Rademacker c , 43 A. Rahaman, 48 V. Ramakrishnan, 60 N. Ranjan, 49 I. Redondo, 32 B. Reisert, 18 V. Rekovic, 38 
P. Ronton, 43 M. Rescigno, 52 S. Richter, 27 F. Rimondi*, 6 L. Ristori, 47 A. Robson, 22 T. Rodrigo, 12 T. Rodriguez, 46 
E. Rogers, 25 S. Rolli, 57 R. Roser, 18 M. Rossi, 55 R. Rossin, 11 P. Roy, 34 A. Ruiz, 12 J. Russ, 13 V. Rusu, 18 
H. Saarikko, 24 A. Safonov, 54 W.K. Sakumoto, 50 O. Salto, 4 L. Santi 1 ", 55 S. Sarkar 11 , 52 L. Sartori, 47 K. Sato, 18 
\. Savoy-Navarro, 45 T. Scheidle, 27 P. Schlabach, 18 A. Schmidt, 27 E.E. Schmidt, 18 M.A. Schmidt, 14 M.P. Schmidt*, 61 
M. Schmitt, 39 T. Schwarz, 8 L. Scodellaro, 12 A.L. Scott, 11 A. Scribano r , 47 F. Scuri, 47 A. Sedov, 49 S. Seidel, 38 
Y. Seiya, 42 A. Semenov, 16 L. Sexton-Kennedy, 18 A. Sfyrla, 21 S.Z. Shalhout, 59 T. Shears, 30 P.F. Shcpard, 48 
D. Sherman, 23 M. Shimojima', 56 S. Shiraishi, 14 M. Shochct, 14 Y. Shon, 60 I. Shreyber, 37 A. Sidoti, 47 P. Sinervo, 34 
A. Sisakyan, 16 A.J. Slaughter, 18 J. Slaunwhite, 40 K. Sliwa, 57 J.R. Smith, 8 F.D. Snider, 18 R. Snihur, 34 A. Soha, 8 
S. Somalwar, 53 V. Sorin, 36 J. Spalding, 18 T. Spreitzer, 34 P. Squillacioti r , 47 M. Stanitzki, 61 R. St. Denis, 22 
B. Stelzer 9 O. Stelzer-Chilton. 43 D. Stentz. 39 J. Stroloe-as 38 D. Stuart,. 11 J.S. Suh 28 A. Sukha.nov 19 



35 

46 



S. Somalwar, 53 V. Sorin, 36 J. Spalding, 18 T. Spreitzer, 34 P. Squillacioti r , 47 M. Stanitzki, 61 R. St. Dc 
B. Stelzer, 9 O. Stelzer-Chilton, 43 D. Stentz, 39 J. Strologas, 38 D. Stuart, 11 J.S. Suh, 28 A. Sukhano^ 
I. Suslov, 16 T. Suzuki, 56 A. TafTard d , 25 R. Takashima, 41 Y. Takeuchi, 56 R. Tanaka, 41 M. Tecchio, 
P.K. Teng, 1 K. Terashi, 51 R.J. Tesarek, 18 J. Thorn^', 18 A.S. Thompson, 22 G.A. Thompson, 25 E. Thomson, 46 
P. Tipton, 61 V. Tiwari, 13 S. Tkaczyk, 18 D. Toback, 54 S. Tokar, 15 K. Tollefson, 36 T. Tomura, 56 D. Tonclli, 18 
S. Torre, 20 D. Torretta, 18 P. Totaro" 1 , 55 S. Tourncur, 45 Y. Tu, 46 N. Turini r , 47 F. Ukegawa, 56 S. Vallecorsa, 21 
N. van RemorteP, 24 A. Varganov, 35 E. Vataga 5 , 47 F. Vazquez^, 19 G. Velev, 18 C. Vellidis, 3 V. Veszpremi, 49 
M. Vidal, 32 R. Vidal, 18 I. Vila, 12 R. Vilar, 12 T. Vine, 31 M. Vogel, 38 I. Volobouev , 29 G. Volpi«, 47 F. Wiirthwcin, 10 

P. Wagner, 2 R.G. Wagner, 2 R.L. Wagner, 18 J. Wagner-Kuhr, 27 W. Wagner, 27 T. Wakisaka, 42 R. Wallny, 9 
S.M. Wang, 1 A. Warburton, 34 D. Waters, 31 M. Weinberger, 54 W.C. Wester III, 18 B. Whitehouse, 57 D. Whiteson^, 4 ' 
A.B. Wicklund, 2 E. Wicklund, 18 G. Williams, 34 H.H. Williams, 46 P. Wilson, 18 B.L. Winer, 40 P. Wittich^, 18 
S. Wolbers, 18 C. Wolfe, 14 T. Wright, 35 X. Wu, 21 S.M. Wynne, 30 A. Yagil, 10 K. Yamamoto, 42 J. Yamaoka, 53 
U.K. YangY 4 Y.C. Yang, 28 W.M. Yao, 29 G.P. Ych, 18 J. Yoh, 18 K. Yorita, 14 T. Yoshida, 42 G.B. Yu, 50 I. Yu, 28 
S.S. Yu, 18 J.C. Yun, 18 L. ZanchV, 52 A. Zanetti, 55 I. Zaw, 23 X. Zhang, 25 Y. Zheng 6 , 9 and S. Zucchelli* 6 

(CDF Collaboration*) 

1 Institute of Physics, Academia Sinica, Taipei, Taiwan 11529, Republic of China 
2 Argonne National Laboratory, Argonne, Illinois 604-39 
3 University of Athens, 157 71 Athens, Greece 
^Institut de Fisica d'Altes Energies, Universitat Autonoma de Barcelona, E-08193, Bellaterra (Barcelona) , Spain 

5 Baylor University, Waco, Texas 76798 
''Istituto Nazionale di Fisica Nucleare Bologna, t University of Bologna, 1-4-0127 Bologna, Italy 
Brandeis University, Waltham, Massachusetts 02254 
University of California, Davis, Davis, California 95616 
9 University of California, Los Angeles, Los Angeles, California 90024 
10 University of California, San Diego, La Jolla, California 92093 
11 University of California, Santa Barbara, Santa Barbara, California 93106 
12 Instituto de Fisica de Cantabria, CSIC- University of Cantabria, 39005 Santander, Spain 
13 Carnegie Mellon University, Pittsburgh, PA 15213 
^Enrico Fermi Institute, University of Chicago, Chicago, Illinois 60637 
15 Comenius University, 842 48 Bratislava, Slovakia; Institute of Experimental Physics, 040 01 Kosice, Slovakia 
16 Joint Institute for Nuclear Research, RU-141980 Dubna, Russia 

11 Duke University, Durham, North Carolina 27708 
18 Fermi National Accelerator Laboratory, Batavia, Illinois 60510 
19 University of Florida, Gainesville, Florida 32611 
20 Laboratori Nazionali di Frascati, Istituto Nazionale di Fisica Nucleare, 1-00044 Frascati, Italy 
21 University of Geneva, CH-1211 Geneva 4, Switzerland 
22 Glasgow University, Glasgow G12 8QQ, United Kingdom 

23 Harvard University, Cambridge, Massachusetts 02138 
2J> Division of High Energy Physics, Department of Physics, 
University of Helsinki and Helsinki Institute of Physics, FIN-00014, Helsinki, Finland 
25 University of Illinois, Urbana, Illinois 61801 
26 The Johns Hopkins University, Baltimore, Maryland 21218 
2 '' ' Institut fur Experimentelle Kernphysik, Universitat Karlsruhe, 76128 Karlsruhe, Germany 



28 Center for High Energy Physics: Kyungpook National University, 
Daegu 702-701, Korea; Seoul National University, Seoul 151-742, 
Korea; Sungkyunkwan University, Suwon 440-746, 
Korea; Korea Institute of Science and Technology Information, Daejeon, 
305-806, Korea; Chonnam National University, Gwangju, 500-757, Korea 
29 Ernest Orlando Lawrence Berkeley National Laboratory, Berkeley, California 94720 
30 University of Liverpool, Liverpool L69 7ZE, United Kingdom 
31 University College London, London WC1E 6BT, United Kingdom 
S2 Centro de Investigaciones Energeticas Medioambientales y Tecnologicas, E-28040 Madrid, Spain 
33 Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 

' 14 Institute of Particle Physics: McGill University, Montreal, 
Canada H3A 2T8; and University of Toronto, Toronto, Canada M5S 1A7 
35 University of Michigan, Ann Arbor, Michigan 48109 
36 Michigan State University, East Lansing, Michigan 48824 
37 Institution for Theoretical and Experimental Physics, ITEP, Moscow 117259, Russia 
38 University of New Mexico, Albuguerque, New Mexico 87131 
39 Northwestern University, Evanston, Illinois 60208 
J, ° The Ohio State University, Columbus, Ohio 43210 
41 Okayama University, Okayama 700-8530, Japan 
42 Osaka City University, Osaka 588, Japan 
43 University of Oxford, Oxford 0X1 3RH, United Kingdom 
''Istituto Nazionale di Fisica Nucleare, Sezione di Padova-Trento, u University of Padova, 1-35131 Padova, Italy 
45 LP N HE, Universite Pierre et Marie Curie/IN2P3-CNRS, UMR7585, Paris, F-75252 France 
46 University of Pennsylvania, Philadelphia, Pennsylvania 19104 
41 Istituto Nazionale di Fisica Nucleare Pisa, 9 University of Pisa, 
r University of Siena and 3 Scuola Normale Superiore, 1-56127 Pisa, Italy 
48 University of Pittsburgh, Pittsburgh, Pennsylvania 15260 
49 'Purdue University, West Lafayette, Indiana 47907 
50 University of Rochester, Rochester, New York 14627 
51 The Rockefeller University, New York, New York 10021 
52 Istituto Nazionale di Fisica Nucleare, Sezione di Roma 1, 
v Sapienza Universita di Roma, 1-00185 Roma, Italy 
53 Rutgers University, Piscataway, New Jersey 08855 
54 Texas A&M University, College Station, Texas 77843 
55 Istituto Nazionale di Fisica Nucleare Trieste/ Udine, w University of Trieste/ Udine, Italy 
56 University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
51 Tufts University, Medford, Massachusetts 02155 

58 Waseda University, Tokyo 169, Japan 
59 Wayne State University, Detroit, Michigan 48201 
60 University of Wisconsin, Madison, Wisconsin 53706 
61 Yale University, New Haven, Connecticut 06520 

We search for b — > s/i + /i~ transitions in B meson (B + , B°, or B®) decays with 924 pb -1 of 
pp collisions at ^fs = 1.96 TeV collected with the CDF II detector at the Fermilab Tevatron. We 
find excesses with significances of 4.5, 2.9, and 2.4 standard deviations in the B + — ► fi + K + , 
B° — > fi + iT K* (892)° , and B° — > fi + (i~(f> decay modes, respectively. Using B — > J/tph (h — K + , 
if* (892) , <f>) decays as normalization channels, we report branching fractions for the previously 
observed B+ and B° decays, B{B+ -> = (0.59 ± 0.15 ± 0.04) x 10~ 6 , and B(B° -> 

li+pTK* (892)°) = (0.81 ± 0.30 ± 0.10) x 10~ 6 , where the first uncertainty is statistical, and the 
second is systematic. These measurements are consistent with the world average results, and are 
competitive with the best available measurements. We set an upper limit on the relative branching 

fraction ^^giP^j^^ < 2.6(2.3) x 10~ 3 at the 95(90)% confidence level, which is the most stringent 
to date. 
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The decay of a b quark into an s quark and two 
muons (b — > sfi + [i~) is a flavor-changing neutral cur- 
rent process, forbidden at tree level in the standard 
model (SM) but allowed through highly suppressed in- 
ternal loops. New physics could manifest itself in a 
larger branching fraction, a modified dimuon mass dis- 
tribution, or angular distributions of the decay products 
different from that predicted by the SM [1-4]. In this 
paper, we report branching ratio measurements of exclu- 
sive decays, where the s quark hadronizes into a single 
meson. Reconstructing a final state with only three or 
four charged final state particles results in smaller back- 
grounds than those expected in a search for inclusive 
B — > X s fi + /j,~ decays. The rare decays, B + — » fi + fi~K + 
and B° ii+^K*(892)° have been observed at the 
B factories [5, 6], with branching fractions of O(10 -6 ), 
consistent with SM predictions [7-14]. The analogous 
decay in the B® system, — > /i + /i~</>, has a predicted 
branching ratio of 1.6 x 10 ~ 6 [15], but has not yet been 
observed [16, 17]. We report branching ratio measure- 
ments from 924 pb" 1 of Collider Detector at Fcrmilab 
(CDF) Run II data that pave the way for future stud- 
ies of larger datasets from which kinematic distributions 
may be measured. We search for B — > n + [i~h decays, 
where B stands for B + , B°, or B®, and h stands for K + , 
K*(892)°, or 0, respectively. The if* (892)°, referred to 
as K*° throughout this paper, is reconstructed in the 
K*° — > K + n~, decay mode, and the cj) meson is recon- 
structed as 4> — > K + K~ . We measure branching ratios 
relative to B — > J/iph decays, followed by J/ip — > 
decays, resulting in the same final state particles as the 
rare decay modes. Many systematic uncertainties cancel 
in the relative branching ratios: 



B(B -> n+fi-h) 
B(B -> J/t(;h) 



x S(J/V> -» 



where N fl + fl - h is the observed number of B — > h de- 

cays, Nj/^ h is the observed number of B — > J/iph decays, 
while ej/^h an d t^+^-h are the combined selection effi- 
ciency and acceptance of the experiment for B — > J j iph 
and B — > n + [i~h respectively. Throughout this report, 
charge conjugate modes are implicitly included. 

CDF II is a general purpose detector, located at the 
Tevatron pp collider [18]. Charged particle trajectories 
(tracks) are detected by the tracking system comprised of 
a seven-layer double-sided silicon microstrip detector and 
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a drift chamber, both in a 1.4 T axial magnetic field. The 
silicon detector [19] ranges in radius from 1.3 to 28 cm, 
and has a single- hit resolution of approximately 15 /im. 
The drift chamber [20] provides up to 96 measurements 
from radii of 40 to 137 cm with a single- hit resolution of 
approximately 180 /im. A muon chamber identification 
system of plastic scintillators and drift chambers [21] is 
located on the exterior of the detector with a central 
part covering \rj\ < 0.6, and an extended region cover- 
ing 0.6 < \r]\ < 1.0, where rj is the pseudorapidity [22]. 
In the central(extended) region, muons are detected if 
their transverse momentum component, px, is greater 
than 1.5(2.0) GeV/c. Events are selected with a three- 
level trigger system. The first trigger level requires the 
presence of two charged particles with px > 1.5GeV/c 
< 0-6) or p T > 2.0GeV/c (0.6 < \rj\ < 1.0), matched 
to track segments in the muon chambers to form muon 
candidates. At the second level, a more restrictive se- 
lection is made by requiring that the muon candidates 
have opposite charge, and that their opening angle in the 
plane transverse to the beamline is less than 120°. At the 
third trigger level, the event is fully reconstructed. The 
trajectories of the muon candidates in the silicon detec- 
tor are required to intersect at a point which is displaced 
transversely from the beamline by at least 100 /zm. 

The offline selection begins with a suppression of ran- 
dom combinations of tracks satisfying the selection re- 
quirements (combinatoric background), by requiring that 
all tracks have px > 0.4GeV/c and match hits from at 
least three layers of the silicon detector. The trajecto- 
ries of a pair of muon candidates that satisfy the trigger 
requirements, and the tracks that form the hadron can- 
didate, are fitted with the constraint that they originate 
from a single vertex in 3-dimensional space to form a 
B candidate. The x 2 probability of the fit is required 
to be greater than 10~ 3 . A K + candidate is a track 
assigned the charged kaon mass, and the K*° (cj)) candi- 
dates are formed from oppositely charged pairs of tracks 
whose invariant mass must lie within 50 (10) MeV/c 2 of 
the K*° (0) mass. For all particles, we use the world av- 
erage values tabulated in Ref. [23] . The ambiguity of the 
mass assignment in B° — ► /j, + /j,~K*° decays is handled 
by choosing the combination whose K + ir~ mass is closer 
to the K*° mass. In reconstructing the B candidates, the 
meson containing a strange quark, h, is required to have 
Pr(h) > l.OGeV/c, and the B candidate is required to 
have pr(B) > 4.0GeV/c. We require that the distance 
of closest approach between the flight path of the B can- 
didate and the beamline, \do(B)\, is less than 120 /jm, 
to reduce the combinatoric background with no loss of 
signal. 

Normalization mode candidates are identified by hav- 
ing a dimuon invariant mass within 50MeV/c 2 of the 
J/tp mass, yielding approximately 12 000 B + — > J /%pK + , 
4400 B° -> J/ipK* , and 800 B° s J/i/><f> decays. The 
B candidate mass distributions of all modes are com- 
patible with a Gaussian of width a = 20MeV/c 2 . To 
reduce backgrounds from B decays to mesons contain- 



ing c quarks, several vetoes are applied to B — > 
candidates, listed below. We eliminate candidates with a 
dimuon mass near the J/V> (,) : 2.9 < m^n < 3.2GeV/c 2 
or 3.6 < m w < 3.75GeV/c 2 respectively. B -> J/ip^h 
decays followed by the radiative decay of the J/ip('> into 
two muons and a photon that is not reconstructed, may 
have a dimuon mass that passes the J/tp^ veto. We re- 
ject these events making use of the correlation between 
the candidate's invariant mass and the dimuon invariant 
mass, and require |(m M(1 A - m B ) - (m w - m^))| > 
100MeV/c 2 . S -> J/tjj^h decays with one hadron 
misidentified as a muon form a potential background to 
the rare decay search. We reject this class of background 
by requiring that all combinations of tracks of the can- 
didate have an invariant mass that differs by at least 
40MeV/c 2 from the J/tpW mass. We also reject candi- 
dates with track pairs with a mass within ±25MeV/c 2 
of the D° — > K~tt + decay, or track triplets compatible 
with a mass within ±25MeV/c 2 of the D+ —> K^tt+tt + 
or — > K + K~w + decays. 

We further improve the signal to background ratio by 
an optimization of the selection based on discriminating 
variables. For this purpose, Monte Carlo (MC) simula- 
tions of the signal processes are used. We generate sin- 
gle b hadrons using the transverse momentum spectrum 
from B — * J/tpX, measured by CDF [18]. Decays of all b 
hadrons are simulated using evtgen [24] , with lifetimes 
from Ref. [23], allowing for a decay width difference be- 
tween the B° s mass eigenstates, AT/T = 0.12 ± 0.06 [25]. 
The decays B° — » J/ipK*° and B® — > J/ip(j) are simu- 
lated according to the polarization amplitudes measured 
by CDF [26]. The dynamics of the rare decay pro- 
cesses are simulated according to the calculations of Ali 
et al. [11]. For B® — > fi + fi~(f), we assume a mixture of 
50% CP-even and 50% CP-odd states. The interactions 
of final state particles are simulated using a GEANT [27] 
model of the CDF II detector, digitized into the CDF 
event format, and reconstructed using the same software 
as in the processing of collision data. The detector sim- 
ulation includes a full emulation of the CDF trigger sys- 
tem. We find that the transverse momentum spectrum of 
the B mesons of this measurement are on average higher 
than in the simulation. This can be explained by the 
presence of a B baryon component with low transverse 
momentum in the measurement used as input spectrum 
of the simulation [18]. To correct for this, the simulated 
events are weighted by a polynomial function, which is 
obtained from the ratio of data to MC transverse mo- 
mentum distribution of the B meson. 

We tighten the candidate selection to obtain the small- 
est expected statistical uncertainty on the measurement 
of the B — > /.i + /j,~h branching ratios by finding the selec- 
tion with the largest value of N sig / N sig + Nbkg ■ The 
expected number of signal events, N sig , is determined 
by scaling the MC yields of the rare decay modes to 
the yields of the normalization channels in the data, 
and correcting by the relative branching ratios and se- 
lection efficiencies. This procedure requires estimates of 



the branching ratios of the rare decay modes. For the 
B + — > /.i + /.i~K + and B° — > fi + /.i~K*° decays we use 
the world average measured branching ratios [23], while 
for B® — > /i + \x~ 4> we use the theoretical estimate [15]. 
The expected number of background events, N^g, is es- 
timated by extrapolating the candidate yield with an 
invariant mass (180 — 300) McV/c 2 higher than the B 
mass, to the signal region. Candidates with an invari- 
ant mass lower than the B mass arc not suitable for 
background estimates, since they contain partially re- 
constructed B decays, not expected in the signal region. 
We find good discriminating power between signal and 
background from the following three quantities: t/a t , a, 
and I. The significance of the proper decay time, t/a t , 
is defined as the proper decay time of the B candidate, 
divided by its uncertainty. The angle a is defined as 
the difference in angle between the B candidate's mo- 
mentum vector and the vector from the primary vertex 
to the [x\xh vertex. The isolation, /, is defined as the 
transverse momentum carried by the B meson candidate 
divided by the transverse momentum of all tracks in a 
cone of AR = yj Arf + A<fi 2 = 1.0 around the direc- 
tion of the B meson candidate, including those of the 
B candidate itself. Here Aij is the difference in pseudo- 
rapidity of the B candidate and each track, and A(j> is 
the difference in their azimuthal angles. Scanning dif- 
ferent combinations of selection thresholds, we find that 
the optimal values are very similar for the three rare 
decay modes: t/at > 14, a < 60mrad, and I > 0.6. 
Applying the optimized selection requirements to the 
normalization channels yields the distributions shown in 
Fig. 1. The shape of the combinatoric background is 
estimated from a sample with poor vertex quality (x 2 
probability < 10~ 3 ). The invariant mass distribution of 
this poor vertex quality sample is fitted with a Gaus- 
sian distribution for the remaining signal contribution 
and an exponential plus a constant to model the back- 
ground. We estimate the background in the signal region 
of the tight selection using this functional form, normal- 
ized to the number of candidates with invariant mass 
(60 — 180) McV/c 2 higher than the B mass. Using this 
method, we find 6361 ± 82 B+ -> J/ifjK+, 2423 ± 52 
B° -> J/iPK*°, and 431 ± 22 B° s -> J/^cj) decays. The 
invariant mass distributions for the rare decay modes are 
shown in Fig. 2. In the ±40McV/c 2 window around the 
B mass we find 90 B+ -> f i + ^-K+, 35 B° -> /i+fi-R* , 
and 11 B® —> /i + fj,~<f> candidates. From simulation stud- 
ies we verified that the present selection accepts decays 
for every kinematically possible dimuon invariant mass, 
with an efficiency difference not exceeding a factor two, 
except for the windows around the J/ip^ that have been 
explicitly excluded. 

For the signal modes, we evaluate three sources of 
background: the combinatoric background, hadrons 
misidentified as muons, and hadrons with misassigned 
mass. As in the normalization mode, the shape of the 
combinatoric background is estimated from a sample 
with poor vertex quality (x 2 probability < 10~ 3 ). The 
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FIG. 1: Invariant mass spectra of (a) J/if)K + , (b) J/ipK*° and (c) J/V></> candidates after applying optimized selection require- 
ments. 



invariant mass distribution of these candidates is fitted to 
an exponential plus a constant. This functional form is 
used to extrapolate the number of candidates with invari- 
ant mass (60 — 180) MeV/c 2 higher than the B mass to 
the signal region. The resulting backgrounds, and their 
statistical uncertainties, consist of 44.3 ± 5.8, 16.3 ± 3.6, 
and 3.1 ± 1.5 candidates for the B + — > /i + /i~A' + , 
B° — > fi + fi~K*°, and B° — » fi + \i~ <f) decay modes, re- 
spectively. Backgrounds from hadrons misidentified as 
muons are estimated from simulation. Final state parti- 
cles from simulated charmless B decays are weighted with 
the muon misidentification probabilities. The misiden- 
tification probabilities for charged kaons and pions are 
measured from a sample of D* + — > D°ir + followed by 
D° — > K~n + decays. We estimate a net background 
of 1.0 events for the B + — > (i + [j,~K + channel, primar- 
ily coming from misidentified B + — > K + Tr~Tt + decays. 
For the B° and B® modes, this background is negligible. 
From simulation we determine the background contribu- 
tion that arises from assigning the wrong masses to the 
hadrons. We estimate 0.4 B° — > *° decays which 

are reconstructed as B° s — > and 0.2 B® — > 

decays reconstructed as B° — > f.i + fi~K*°. 

We calculate the signal yield by subtracting the pre- 
dicted background from the number of candidates in the 
signal window. We find an excess in the B signal region 
in all three channels, and determine the significance by 
calculating the Poisson probability for the background 
to fluctuate to the number of observed events or higher, 
taking into account the uncertainty on the background. 
We find an equivalent Gaussian significance of 4.5, 2.9, 
and 2.4 standard deviations, respectively for the B + , B°, 
and B° modes. 

Table I lists the systematic uncertainties associated 
with the relative selection efficiency, discussed in fur- 
ther detail below. To take into account the uncertainty 
on the relative efficiency due to uncertainties in the dy- 
namics of the rare decays, we apply the weak form fac- 
tors from Ref. [28, 29] and evaluate the largest differ- 
ence between the reconstruction efficiency and the cen- 



tral value. We find an uncertainty of 3.1% or less. The 
uncertainty related to the pt(B) spectrum is evaluated 
from the change in relative efficiency when the three pa- 
rameters in the pt{B) weighting function are varied by 
one standard deviation of their values determined from 
fits to data, taking into account correlations. The dif- 
ference reaches 1.4% in the B® mode. The muon trigger 
efficiency close to the 1.5GeV/c pr threshold is poorly 
known. We find a change in the relative efficiency reach- 
ing 1.3% if the minimum pr is varied by ±100MeV/c, 
covering the range over which the trigger efficiency rises 
from zero to one. The final state particles of the rare 
decay modes have approximately 10% more low momen- 
tum tracks than those of the normalization channels. The 
simulation models the track reconstruction efficiency to 
an accuracy of 2% in the pr range of 0.4 — 1.5 GeV/c. We 
therefore assign a systematic uncertainty on the relative 
efficiency of 2%xl0%=0.2%. We find an uncertainty of 
8.3% on the relative efficiency of i?° decays due to the 
unknown fraction of the short-lived CP-even state in the 
rare decay mode. We assume that the CP-even frac- 
tion is 0.5, and assess systematic uncertainties for the 
extremes of the fraction valued at 1.0 and 0.0. In addi- 
tion, the uncertainty on AT/T contributes another 2.6%, 
resulting in a total uncertainty of 8.7% associated with 
the B® decay width difference. This is the largest sys- 
tematic uncertainty on the relative efficiency of the B® 
mode. We evaluate the effect of the uncertainty on the 
fraction of J / tp mesons produced with a longitudinal po- 
larization by varying the fraction measured at CDF [26] 
by ±lcr. The effect is 0.6% for B° -> J/tpK* and 0.1% 
for B° s -> J/ip<f>. 

The statistical uncertainties of yields in the normal- 
ization channels are included as systematic uncertainties, 
which range from 1.3% for the B + channel to 5.1% for the 
B® channel. Cabibbo-suppressed B + — > J/tfjn + decays 
contribute 0.1% to the yield of the normalization channel 
for the B + channel. We introduce the full effect into our 
estimate of the systematic uncertainty without correct- 
ing the result. The relative efficiencies between signal 



TABLE I: Systematic uncertainties on the relative efficiency 
quoted in percent. 



Channel B + B u ~BT 

Theory model 1.5 3.1 1.6 

p T (B) spectrum 0.6 1.3 1.4 

Trigger turn-on 1.3 1.3 1.2 

Low momentum hadrons 0.2 0.2 0.2 

5° decay width difference — — 8.7 

Polarization — 0.6 0.1 

Norm, channel statistics 1.3 2.1 5.1 

B + — > J/ipir + contribution 0.1 — — 

MC statistics 1.6 2.6 2.2 

Total 2.9 5.0 10.6 



TABLE II: Summary of systematic uncertainties quoted in 
percent. 



Channel 


B+ 


B u 


B s 


Total rel. eff. uncertainty 


2.9 


5.0 


10.6 




1.0 


1.0 


1.0 


Background prediction 


5.2 


3.1 


10.2 


B(B -> J/i>h) 


3.5 


4.5 


35.5 



and normalization channels, of 0.71 ± 0.01, 0.74 ± 0.02, 
and 0.84 ± 0.02 for the B+ , B°, and B° decays, respec- 
tively, have uncertainties reaching 2.6%, that arise from 
the finite size of the MC samples. 

The predicted background values depend on the shape 
of the background function. We evaluate the change in 
the background yield calculation when using a sample 
that is similar to the one resulting from the optimal se- 
lection, but has less stringent requirements on the three 
optimization variables instead of the default sample with 
poor vertex quality. We compare two functional models 
for the background description, the default exponential 
plus a constant, and a simpler linear extrapolation func- 
tion. We find a systematic uncertainty on the background 
prediction reaching 10% for the B® channel. 

To calculate the systematic uncertainty on the rela- 
tive branching fraction, the relative efficiency uncertainty 
is summed in quadrature with the uncertainty of the 
J '/ip — > branching ratio, and with the systematic 

uncertainty on the number of signal candidates. 

To calculate the absolute branching fractions, we use 
the world average branching fractions of the normaliza- 
tion channels [23]. These branching fractions have un- 
certainties of 3.5%, 4.5%, and 35.5%, respectively, for 
B + , B°, and B®, which are added in quadrature to the 
systematic uncertainties on the relative branching ratios. 
A summary of the systematic uncertainties is given in 
Table II. 

We calculate relative branching ratios using Eq. (1), 
and find B(B+ -> n+ yT K+)/B(B+ -> J/4>K+) = 
(0.59±0.15±0.03) x 10~ 3 , B(B° fj+ fj,~ K* )/ 'B(B° -> 



J/^K*°) = (0.61 ± 0.23 ± 0.07) x 10~ 3 , and B{B° -> 
fi+fj,-<t>)/B(B° -> J/^4>) = (1.23 ± 0.60 ± 0.14) x 10~ 3 , 
where the first uncertainty is statistical, and the sec- 
ond is systematic. We use the world-average branch- 
ing ratios of the B° and B + normalization channels [23] , 
resulting in the following absolute branching fractions: 
B(B+ -> H+LI-K+) = (0.59 ±0.15 ±0.04) x 10~ 6 , and 
B(B° -> n+fi-K* Q ) = (0.81 ± 0.30 ± 0.10) x 10~ 6 . To 
obtain an absolute branching ratio for the B® — > fi + fj,~<fi 
decay, we use B{B° S -> J/ip(p) = (1.38 ± 0.49) x 10~ 3 , 
obtained from correcting the CDF measurement [30] for 
the current value of f s /fd [23], the B° s to B° produc- 
tion ratio, resulting in B(B° S -> fi + ^(f>) = (1.70 ± 0.82 ± 
0.64) x 10~ 6 . We find a significance of only 2.4 standard 
deviations for the B® — > /i + /i~(/> decay mode. There- 
fore, we choose to set a limit on this decay. We use a 
Bayesian integration assuming a flat prior [31], and find 
B(B° S -» n+fi-(f))/B(B s -» < 2.6(2.3) x 10~ 3 at 

the 95(90)% confidence level (C.L.). We also set an upper 
limit on the absolute branching ratio, taking into account 
the uncertainty of the branching ratio of the normaliza- 
tion channel [32]. We assume that the prior probability 
density function representing the uncertainty in the nor- 
malization channel is a log normal distribution with a 
mean equal to the central value, and a width parame- 
ter equal to the quoted uncertainty of 35.5%. We obtain 
B(B° -> Ai+M~0) < 6.0(5.0) x 10~ 6 at 95(90)% C.L. The 
main ingredients and the results are summarized in Ta- 
ble III. 

In conclusion, our measurements of the B + and B° rare 
decay modes are consistent with the SM predictions, and 
with previous measurements [5, 6]. The relative limit on 
B® — > f-i + n~(f) is consistent with the SM predictions, and 
is the most stringent to date. 
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FIG. 2: Invariant mass spectra of (a) fi + K + , (b) fj,~ K*° and (c) /i + fj,~4> candidates. The superimposed curves are a sum 
of a single Gaussian with a width of 20 MeV/c 2 representing the signal, and a curve representing the background as determined 
by the procedure described in the text. The curves are not drawn for masses below the signal window, since they are not 
expected to predict the background level where partially reconstructed B decays contribute. 



TABLE III: Summary of the main ingredients and results of this analysis. Where one uncertainty is quoted, the uncertainty is 
of statistical nature. Where two uncertainties are quoted, the first is statistical, and the second is systematic. 

"Decay mode B + -> g ^ K + B u -> fi + ^ A'* u fi + n~ (j> 

IV^ 90 35 11 

N bkg 45.3 ±5.8 16.5 ±3.6 3.5 ±1.5 

N sig 44.7 ±5.8 18.5 ±3.6 7.5 ±1.5 

Gaussian significance 4.5 a 2.9 a 2.4 a 

N jm 6361 ±82 2423 ± 52 431 ± 22 

e M + M -h/ej/*h 0.71 ±0.01 0.74 ± 0.02 0.84 ± 0.02 

Rel B x 10 3 0.59 ± 0.15 ± 0.03 0.61 ± 0.23 ± 0.07 1.23 ± 0.60 ± 0.14 

Abs B x 10 6 0.59 ± 0.15 ± 0.04 0.81 ± 0.30 ± 0.10 1.70 ± 0.82 ± 0.64 

Rel B 95(90)%C.L. limit xlO 3 2.6(2.3) 

Abs B 95(90)%C.L. limit xlO 6 — — 6.0(5.0) 
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